Monochromatic Neutrino Signals from Dark Matter Annihilation 
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Dark matter (DM) annihilations in the Sun to neutrino-antineutrino pairs have striking signatures 
in neutrino detectors such as IceCube and KM3. We make a model independent study of the signals 
after propagation of the neutrinos from the center of the Sun to the Earth. A large spin-dependent 
DM capture cross section in the Sun make the discovery prospects robust. 



INTRODUCTION 

Dark matter is perhaps the best evidence for physics 
beyond the standard model (SM). The amount of dark 
matter is known to be about 22% of the total matter- 
energy budget of the universe, while visible matter only 
comprises of about 4% of the universe Any realistic 
model of physics beyond the SM must supply a DM can- 
didate that can reproduce the observed relic density and 
many models have been suggested that can yield a DM 
candidate. Generally, a parity in the model leads to a 
stable DM particle. The leading DM candidate is a neu- 
tral weakly interacting stable particle that was produced 
in the Big Bang. 

Vigorous experimental efforts are underway to iden- 
tify the nature of DM. A direct way to identify the DM 
particle is to observe the recoil of a nuclear target. A 
spin-independent (SI) interaction is probed by detecting 
the recoil of DM from high mass nucleus; the coherent 
interaction increases the total cross section substantially 
compared with light nuclei. The best current limit comes 
from the XENON10 experiment which found an upper 
bound of 4.5 x 1(T 8 pb for a DM mass of w 30 GeV [|; it 
is expected that CDMS will soon improve on their current 
limit of 1.6 x 1(T 7 pb 3]. The spin-dependent (SD) inter- 
action is governed by the spin of the nucleon. The upper 
limits on the SD cross section are 6 orders of magnitude 
weaker than SI with the best bound from ZEPLIN-II at 
0.07 pb [J]. For a recent discussions of various direct and 
indirect DM searches, see Ref. [1, @, 0) ■ 

At the Large Hadron Collider (LHC), a characteristic 
signature of the DM particle will be missing transverse 
energy. Typically, the DM particle is created at the end 
of a cascade decay chain of new physics particles. The 
reconstruction of the mass of the DM particle at colliders 
is challenging but feasible. 

Additionally, DM is being sought in astrophysics ex- 
periments through annihilations to 7-rays H ToL 
[l2T | , antideuterons [13[ and positrons [3, UH, [l 
galactic halo. 



The method of DM detection on which we concentrate 
in this study is the detection of muon neutrinos from 
DM annihilations in the Sun and Earth with km 2 size 
neutrino detectors [It], EH- Currently, the IceCube ex- 
periment at the South Pole is underway; it has a muon 
energy threshold of about 50 GeV [19|, |20( . The planned 
km 2 size detector in the Mediterranean Sea, known as 



KM3 [2lJ, should detect neutrino interactions above a 10 
GeV muon energy threshold. 

DM is expected to be gravitationally captured in the 
core of the Sun and Earth and subsequently annihilate. 
The rate of capture is strongly dependent on the SD and 
SI scattering rates that are being limited by direct de- 
tection experiments. For DM capture in the Sun, the SD 
cross section is the dominant interaction since the Sun is 
primarily composed of Hydrogen which has a net spin. 
The situation for the Earth is very different as the con- 
tributions to the capture rate via the SD interaction are 
negligible for heavy nuclei. Overall, since the SD cross 
section is nearly 6 orders of magnitude less constrained 
than the SI cross section, models that predict large SD 
rates have greater potential to yield an observable neu- 
trino flux from DM annihilations in the Sun. 

The study of the properties of DM annihilations to neu- 
trino pairs has advantages over a more general neutrino 
spectrum. The injected neutrino spectra from this sim- 
ple process consists of a line at E v = Mn m 1 ■ Models in 
which a neutrino line spectrum from DM annihilations 
can be realized include a Dirac fermion, Kaluza-Klein 
spin-1 boson and a scalar that couples to a new vector 
field. 

The solar capture rate of dark matter in the galactic 
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1 Using this injection spectrum, which can be taken as a delta func- 
tion, one can reproduce any general spectrum after propagation 
via the Green's function of the evolution equation. 
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halo is approximately given by [2J 
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where pi OC ai and w; OC a; are the local density and velocity 
of relic dark matter, respectively. The average density of 
Piocai ~ 0.3 GeV/cm 3 may be enhanced due to caustics 
in the galactic plane [2j|. The effective strength of the 
capture cross section of DM with solar matter is given 
by a cap = erf D + <rg, + 0.07crf/. The factor of 0.07 
before agf comes from the relative abundance of helium 
and hydrogen in the Sun, as well as other dynamical and 
form factor effects 24| . The cross sections determine how 
efficiently the Sun slows and captures DM. The value of 
a ca p has considerable variation with models. 

The capture rate of dark matter in the galactic halo 
by the Earth is approximately given by |25| 
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is the reduced mass of the 



1 Mdm?TIv 

where p p - DM = Mdm+ ^ p 
DM-nucleon system and q cap = crgj. The form factor, 
/®(MrjAf), is given in Ref. ,25]. A positive signal from 
the Earth, but not the Sun would strongly suggest that 
the DM particle has no significant SD interaction with 
matter. 



Here m is the neutrino mass matrix in the gauge eigen- 
state basis, E v is the neutrino energy, Gf — 1-66 x 
10~ 5 GeV~ 2 is the Fermi constant and N e (r) and N„{r) 
are the electron and neutron densities in the Sun [2q |. 
The Neutral Current (NC) and Charged Current (CC) 
source terms describe the absorption and re-injection of 
neutrinos caused by NC and CC processes while the in- 
jection term describes the spectra injected from DM an- 
nihilation in the core of the Sun or Earth 



MUON RATE AT ICECUBE 

We calculate the detection rates with up-going muons 
following the simulation for IceCube outlined in Ref. [28| . 
The neutrino flux at the surface of the Earth is given by 
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BF(DMDM -> vv), (5) 



where the parameter C is the DM capture rate in the 
Sun (Co) or Earth (C®) given explicitly above in Eqs. 
((TJ and {2}, respectively. The factor of \ is associated 
with the fact that two DM particles produce one anni- 
hilation event; here R — 1.49 x 10 11 m is the Earth-Sun 
distance for annihilation in the Sun, or R = 6.4 x 10 6 m 
is the radius of the Earth for annihilation in the Earth's 
core. The first factor in Eq. ([5]) describes the flux at 
the surface of the Earth while the j^f- term is the nor- 
malized differential rate and BF (DM DM — > vv) is the 
fraction of annihilations to neutrino-antineutrino pairs. 

The muon rate in a km 2 area detector such as IceCube 
can be determined by folding the neutrino flux with the 
muon production cross section 3, 24 1 



NEUTRINO PRODUCTION AND 
PROPAGATION 

We assume that DM annihilates at the capture rate 
to neutrino pairs with a democratic ratio of v e : v^ : 
v r : v e : : v T = 1 : 1 : 1 : 1 : 1 : 1, producing 
a characteristic line shape in the initial neutrino energy 
spectra. The neutrinos are then propagated through the 
Sun as in Ref. [l8| . The evolution equation is given by Q 
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where p(r, E) is the complex density matrix in the gauge 
eigenbasis describing the state of the neutrino of energy 
E at a distance r from the center of the Sun. The Hamil- 
tonian, H , includes the effects of vacuum oscillation from 
nonzero neutrino mass splitting and the matter interac- 
tion: 
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The densities of protons and neutrons near the detec- 
tor are p v = \Na cm -3 and p n = ^Na cm" 3 , respec- 
tively, where is Avagadro's number 2 . The muon 
range, R^E^), is the distance a muon travels before its 
energy falls below a threshold energy, Ej^ r . We take 
50 GeV, which is optimistic for IceCube, but the 



j^thr 

muon thresholds are expected greatly improve 2jJ. Due 
to the long muon range, the fiducial volume of the de- 
tector can be factored as the product of the range and 
the cross sectional area of the detector, called the effec- 
tive area. The effective area of the detector, A e ff, is 
calculated for IceCube following Ref. [28j |. 



2 Since the muon range is at most 1 km for a 1 TeV muon, the 
point of muon production can be assumed to be in ice, rather 
than the Earth's crust. 
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We apply Eq. (J6]) to calculate neutrino signals in Ice- 
Cube from DM annihilations and the backgrounds from 
atmospheric neutrinos. The backgrounds from atmo- 
spheric neutrinos are calculated following Ref. 

Annihilations in the Sun: The DM accumulation and 
annihilation rates in the Sun are assumed to equilibrate: 
r anli = \Cq. In Fig. QJi, we show the signal muon 
event rate and atmospheric background in IceCube 
for four DM masses. The numbers of predicted signal 
events, based on a SD capture cross section of a cap = 
10~ 3 pb and a branching fraction of 10%, for Mdm values 
of 100, 200, 400 and 1000 GeV are 1800, 3000, 11000 
and 160 events per year, respectively. The atmospheric 
background for the 50 - 300 GeV window is 10.4 events 
per year. 

We show the statistical significances in Fig. [TJd for 
a broad energy region of E l ^ r < < 300 GeV. This 
choice of the upper cut at 300 GeV is because the large 
attenuation of high energy neutrinos by the CC and NC 
interactions in the Sun suppresses signal contributions at 
higher energies. Once evidence is seen for a signal, the 
window can be narrowed to improve the significance. 

Given the restrictive experimental limit on o~si , there is 
no hope to see a corresponding signal from annihilation to 
neutrino pairs in the Sun. The SD cross section must be 
large for IceCube to probe the ^-line signal. The Super- 
Kamiokande experimental limit on the induced muon flux 
is approximately 5 x 10~ 15 cm~ 2 s -1 [3(J, and is given 
by the vertical line in Fig. [TJd. 

Information on the nature of the DM can be extracted 
from the muon spectra of the signal shown in Fig. [1] 
A determination of the DM mass can be made from the 
muon energy distribution and the capture rate of DM in 
the Sun can be determined from the total signal rate. 
The capture rate in turn gives information on the local 
galatic DM density and its velocity (cf. Eq. ([1])) that 
can complement inferences from direct detection experi- 
ments. Additional properties of DM can be extracted by 
measuring the neutrino flavor ratios (3~0 ] . 

Shown in Fig. [TJd is the statistical significance (with 
one year of data) of a signal relative to the atmospheric 
neutrino background vs a cap x BF (DM DM -> vv) for 
various DM masses. The 5a discovery regions for Ice- 
Cube are above the blue horizontal lines. The limit 
from XenonlO of 4.5 x 10~ 8 pb conservatively assumes 
BF(DMDM — > vv) = 1 and that the capture rate is 
entirely given by the SI interaction (as in the case of 
scalar DM with no new gauge fields). Therefore, a model 
which has no SD interaction is excluded to the right of 
the XenonlO line. For values of BF(DMDM -> vv) < 1, 
the vertical exclusions lines move to the left. It is evident 
that the prospects for IceCube to discover DM annihila- 
tions to neutrinos from the solar core are robust for the 
SD case, but are dim for the SI case. 

Annihilations in the Earth: We optimistically assume 
^C(9 for which the equilibration time is 
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FIG. 1: (a) Predicted muon event rates in IceCube 
from monochromatic neutrinos. The muon rates assume 
BF(DMDM —> vv) = 0.1 and a capture cross-section a cap 
1 pb that is dominated by the SD interaction. The signal 
region includes < 300 GeV. (b) Statistical significance for 
discovery. The discovery region is given above the blue hor- 
izontal lines at 5er for one and five years of data. The limits 
from XenonlO conservatively assume the maximum value of 
BF(DMDM — » vv) = 1 while the Super-Kamiokande exclu- 
sion is a direct measurement limit on the high energy muon- 
neutrino flux. 



longer than the age of the solar system. Otherwise, 
T = ±C© tanh 2 (y/C^A^t Earth) , where A e = is 
the annihilation rate times relative velocity per unit vol- 
ume and tEarth — 4.5 Gy. However, equilibrium may 
not have yet been reached in the case of the Earth 3 . If 
not in equilibrium, the annihilation rate is dependent on 
the DM annihilation cross section. By assuming the DM 
capture and annihilation processes arc in equilibrium, we 



3 The equilibration time for the Earth is expected to be approxi- 
mately a thousand times longer than that of the Sun [251 . 



are assuming the best-case scenario. 

The angular area of acceptance is larger than the solar 
case due to the less concentrated density profile of the 
dark matter in the Earth. It is estimated [25J that about 
98% of the DM is contained within a window of 
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about 9 Z = 180°. Thus, the statistical significance of the 
annihilation signal improves with larger Mdm- 

In Fig. [2k, we show the expected muon rates in Ice- 
Cube for a cap = 10~ 8 pb and BF (DM DM -> vv) = 
0.1. Constraints from Super Kamiokande, with <!>„ J$ 
3 x 10~ 14 cm~ 2 s _1 [30j |. are not strong enough to limit 
the possibility from a signal from the Earth which we pre- 
dicted to be about $„ 10~ 19 cm~ 2 s _1 . The signal is 
substantially greater than the background; the numbers 
of predicted signal events based on a SI capture cross sec- 
tion for M DM values of 100, 200, 400 and 1000 GeV are 
1.4, 2.6, 1.8 and 0.4 events per year, respectively. The 
signal is enhanced at Mdm ~ 60 GeV due to the abun- 
dance of Iron in the Earth [25| . The high IceCube muon 
energy threshold of 50 GeV limits the opportunity of dis- 
covery. For the KM3 detector 21 1, with a lower energy 
threshold, the prospects for discovery may be more likely. 

In summary, the search for DM annihilations in the 
Sun to neutrino pairs is very promising, provided that 
the SD interaction dominates the DM capture rate. The 
case for corresponding discovery of DM annihilations in 
the Earth is less promising due to the strong limits on 
the capture rate because of the stringent upper bound on 
the SI cross-section. 
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illustrated with a SI capture cross-section of a cap = 10 -8 pb 
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